Low-energy spectral features of supernova (anti) neutrinos in inverted hierarchy 
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In the dense supernova core, self-interactions may align the flavor polarization vectors of v and 
77, and induce collective flavor transformations. Different alignment ansatzes are known to describe 
approximately the phenomena of synchronized or bipolar oscillations, and the split of v energy 
spectra. We discuss another phenomenon observed in some numerical experiments in inverted 
hierarchy, showing features akin to a low-energy split of 77 spectra. The phenomenon appears to 
be approximately described by another alignment ansatz which, in the considered scenario, reduces 
the (nonadiabatic) dynamics of all energy modes to only two v plus two V modes. The associated 
spectral features, however, appear to be fragile when passing from single- to multi-angle simulations. 
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I. INTRODUCTION 

Supernova (SN) neutrinos continue to be a subject of 
great interest in astroparticle physics In particular, 
renewed attention is being paid to collective features of 
flavor transformations induced by v (77) self-interactions 
(see and refs. thereinb The observed collective phe- 
nomena of synchronized (U and bipolar oscillations, 
and the split of v spectra @, @| can be largely under- 
stood, in flavor space, in terms of various alignment (or 
antialignment) ansatzes for the Bloch polarization vec- 
tors of v (P) and 77 (P) at different energy E. Conversely, 
lack of alignment indicates flavor decoherence Q . 

In this work we deal with a minor — yet interesting — 
phenomenon previously observed in the numerical exper- 
iments of Q and then of [H| , akin to a low-energy 
"antineutrino spectral split" in inverted hierarchy. Dif- 
ferently from the neutrino case, this feature appears to 
have a nonadiabatic origin. We show that the 77 spectral 
split observed in Q can be approximately described in 
terms of P and P alignment along four global modes (at 
low and high energy, for v and 77). We also discuss a 
scenario with different input spectra, where the 77 split 
is slightly more evident, at least in the approximation of 
averaged trajectories. The effects described herein might 
play a role, in principle, in future low-energy SN neutrino 
observations via 77 e + p — > n + e + scattering. 

II. EQUATIONS AND NUMERICAL RESULTS 

As in Q, we consider an effective two-flavor {v e ,v x ) 
scenario with mass-mixing parameters (Am 2 , #13). After 
trajectory averaging (single-angle approximation), the P 
and P modes obey the equations of motion (EOM) 

P = H x P = (+uB + Az + fxD) x P , (1) 

P = H x P = (-wB + Az + fxD) x P , (2) 

where ui = Am 2 /2E is the vacuum oscillation frequency, 
A = \/2GfN b is the matter potential related to the e~ 
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FIG. I: Radial profiles adopted for the matter (A) and self- 
interaction (fj,) potentials, in the range r G [10, 250] km. 



density (N e ), /i = ^/2Gf{N+N) is the neutrino-neutrino 
potential related to the effective v and 77 densities (N and 
N), and B = (sin 2#i3, 0, =F cos 2^13), where the upper 
(lower) sign refers to normal (inverted) hierarchy. Only 
the latter case will be considered hereafter. 

Global polarization vectors are defined as integrals 
weighted by v and 77 spectral densities (n and n) 0] 
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and their difference provides the vector D = J — J. 
Numerically, we take Am 2 = 2 x I0" 3 cV 2 , so that 
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We also assume sin #13 = s 2 3 = 10~ 4 as reference value. 
It is then B ~ z, with D z = D-z ~ D-B being a constant 
of motion Q. Concerning the potentials A and /i, their 
radial profile above the neutrinosphcrc (r > 10 km) are 
taken as in Q and shown in Fig. 1. 
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FIG. 2: Single-angle simulation in inverted hierarchy: Final 
fluxes (at r = 250 km, in arbitrary units) for different neutrino 
species as a function of energy. Initial shapes (r = 10 km) are 
shown as dotted lines to guide the eye, with average energies 
reported on top. 



Figure 2 shows typical results at the end of collective 
flavor transformations in inverted hierarchy (not much 
happens in normal hierarchy), as obtained by solving 
Eqs. ([l][2]) with the same initial conditions as in Q. The 
most prominent feature in Fig. 2 is the v spectral split 
in the left panel at a critical energy E c ~ 7 MeV, but a 
"minor" V spectral split is also visible in the right panel 
of Fig. 2 at E c , as observed in @. A similar mi- 
nor feature has been later observed as a small low-energy 
"shoulder" 

[13 or "bump" [13 in the final v spectra. 



III. A DESCRIPTION WITH FOUR VECTORS 

In general, numerical solutions to Eqs. (HJ0) can be 
partly understood via appropriate "alignment" ansatzes 
for P and P modes. In particular, for large /x, one can 
posit that they are separately pinned to their sum (i.e., 
P || J and P || J) 0. The EOM for J and J can then be 
cast in a closed form, formally equivalent to the EOM of 
a gyroscopic pendulum [f| , which elegantly explains the 
so-called synchronized and bipolar oscillations observed 
in numerical experiments [3, [B[ . Alternatively, for slowly 
varying (adiabatic) fi, one can posit alignment with the 
Hamiltonians in Eqs. (|1I2[) : P || H and P |j H, explaining 
nicely [|[ the so-called v spectral split 0, 

The above two ansatzes are, in part, exclusive. The 
"pendulum" solution describes well the synchronized and 
bipolar oscillation phases, but not the final v spectral 
split. Conversely, while the adiabatic solution describes 
well overall flavor evolution features, but not the tran- 
sient, nonadiabatic bipolar oscillations. However, nei- 
ther describe the low-energy 77 spectral split of Fig. 2. 
This phenomenon must thus be related to the breaking 
of global alignment and of adiabaticity. 

Adiabaticity breaking is expected in our scenario, es- 
ecially at low E. According to the criterion proposed in 
(which is independent of both #13 and A), the evo- 
lution becomes nonadiabatic when the value of fi ap- 
proaches a typical vacuum frequency oj. For E < 40 MeV 
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FIG. 3: Single-angle simulation in inverted hierarchy: Moduli 
and z-components of the four vectors defined in Eqs. (|6I7|I . 



(Fig. 2) , the nonadiabaticity criterion is roughly \i ~ cu > 
0.13 km -1 , and is satisfied at any radius r in Fig. 1. 
Moreover, since \x is a decreasing function of r, adia- 
baticity is violated more at high lo (low E). However, 
we have not been able to sharpen this qualitative expec- 
tation and to "predict," for instance, the observed value 
of the v critical energy E c in terms of a specific transi- 
tion to nonadiabaticity. For instance, we find that E c 
depends also on 6i 3 and A; e.g., it decreases by ~2 MeV 
by switching off the matter term A, but is restored to 
E c ~ 4 MeV for sf 3 = 10~ 6 at A = 0. Numerical ob- 
servations thus suggest that the origin of the E c value is 
linked to the complex dynamics close to the z axis, which 
can be rather subtle. 

Here we adopt a pragmatic approach, and take E c and 
E c from the numerical experiment at face value. Our 
main observation is that these energies naturally sepa- 
rate "low" (I) and "high" (h) energy modes, and that 
these modes appear to be separately pinned to their 
global sum. More precisely, we define four new global 
vectors by splitting the energy intervals in Eqs. (|3|4| as 
[0, E c ] U [E c , 00] for v and [0,E C ] U [E c , 00] for V, 



J 
J 



(6) 
(7) 



We find that the moduli of the above vectors (J/ = |J;|, 
etc.) are approximately conserved during the evolution. 

Figure 3 shows, for each of the four vectors defined 
above, both the modulus and the z-componcnt as a func- 
tion of r. The initial moduli, calculated for E c ~ 7.6 MeV 
and E c ~ 3.9 MeV, are: 



Ji -0.149 
Ji ~ 0.008 



J h ~ 0.096 , 
J h ~ 0.101 . 



(8) 
(9) 
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Such values do not change much with r, despite the fact 
that the four vectors, initially aligned in the synchro- 
nized phase (r < 70 km), spread apart and nutate dur- 
ing the bipolar phase (70 < r < 120 km) and finally 
split up with low-energy (high-energy) modes aligned 
(antialigned) with the z axis, corresponding to no (full) 
flavor conversion. 

A posteriori, this behavior is equivalent to posit the 
(approximate) alignment ansatz: P || J; for E < E c , 
P || 3^ for E > E ci and similarly for antineutrinos. Con- 
sistency of this ansatz with conservation of D z from the 
initial (aligned) state to the final (split) state requires 

{■h + Jh) - (Ji + Jh) = (Jj - A) - (Ji - J h ) , (10) 



namely, 



Jh = Jh 



(11) 



which, in our case, is fulfilled within a few % [lJl- The 
fact that the flavor evolution in Fig. 3 is reasonably de- 
scribed by such an ansatz is a clear signal of deviations 
from the purely adiabatic solution, which would instead 
predict an alignment of the P's and P's to linear combi- 
nations of B and D [(|[l3j|, which we do not observe. On 
the other hand, the fact that E c may change somewhat 
with r in the adiabatic solution [l3T | (while we have as- 
sumed r-independent E c and E c values a priori) suggests 
that our ansatz, despite its effectiveness, may be further 
refined. 

Summarizing, the global flavor evolution in our sce- 
nario appears to be effectively described in terms of four 
vectors with nearly conserved lengths, collecting v energy 
modes below and above a critical value E c , and 77 modes 
below and above a (lower) critical E c . It is then useful 
to check if the behavior in Fig. 3 is also captured by the 
solutions of EOM reduced to exactly four modes, namely, 

j; = (+cj/B + Az + fjB) x J, , (12) 

J h = (+oj h B + Az + fxD) x J h , (13) 

j, = (-WiB + Az + /zD) x J ; , (14) 

J h = (-ZUftB + Az + fiD) x J h . (15) 

where the four frequencies are defined as averages of cu 
(with weights n e — n x and n e — n x for v and V, respec- 
tively 0) in the corresponding energy ranges ([0, i? c ] for 
u>i, etc.). Numerically, in our case: u>i ~ 1.21, u>h — 0.69, 
ZJi ~ 2.04, and Zu h ~ 0.64 (all in km" 1 ). The initial 
conditions correspond to alignment of all vectors to +z, 
with moduli given by Eqs. (|8I9[) . Such moduli are triv- 
ially conserved by construction. The z-componcnts are 
more interesting, and evolve in a way rather similar to 
those in Fig. 3 (not shown). In particular, since conser- 
vation of D z prevails over exact alignment, the vectors 
J; and J; (Jh and Jh) are almost — but not exactly — 
aligned (antialigned) to +z. Therefore, even this simpli- 
fied (4-mode) case shows an "antineutrino spectral split" 
to a good approximation, via the opposite behavior of J; 
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FIG. 4: MSW effect on the z-component of v modes at very 
low energies (E = 0.012, 0.154, and 0.444 MeV). 



and Jh, similar to the results of the complete (100-modc) 
spectral case shown in Figs. 2 and 3. 

We conclude by discussing the tiny effect of ordinary 
Mikheev-Smirnov-Wolfenstein (MSW) nonadiabatic res- 
onances in matter for very low energy V modes. In our 
scenario (sf 3 = 10 -4 and [iD z <C A) the V MSW reso- 
nance condition H • z = (i.e., £jcos2#i3 = A + /j,D z ) 
approximately reads 



LO ~ A 



(16) 



This condition is met before the end of collective ef- 
fects (r < 250 km) for uj > 5 km" 1 (E < 1 MeV), 
see Fig. 1. When passing from Aj/w 3> 1 (initial state 
i) to Xf/uj <C 1 (final state /) through the resonance, 
the v e flavor survives with probability P ee ~ P c , where 
P c = exp(— 2iru>Si 3 X/ X) is the level crossing probability 
fl6| . The relation between and the z-components of 
the V polarization vectors P [9( provides the depth of the 
MSW resonance as 



p f z ~pi(2P c -i) 



(17) 



Figure 4 shows the radial profile of P z , as obtained in our 
numerical experiment for three representative (very low) 
energies. The radius and depth of the MSW resonance 
agree with Eqs. (fit!]) and (fTT)) . respectively. The lower the 
energy, the deeper the resonance, the more inverted is the 
final polarization vector P. However, the non-collective 
(uj > /it) MSW effect is limited to such low energies to be 
practically unobservable in Fig. 2, and is numerically ir- 
relevant for the 4-mode approximation described above. 
Finally, we mention that, in normal hierarchy, we find 
that the spectra of v and V are basically unaffected by 
self-interaction effects, the only tiny change being a spec- 
tral swap of neutrinos for E < 1 MeV, induced by MSW 
effects (not shown). 
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FIG. 5: As in Fig. 2 (single angle), but for (E x ) = 18 MeV. 



IV. SINGLE- VERSUS MULTI-ANGLE 

So far we have examined in detail the low E features 
of the same scenario of Q in single-angle approximation. 
It was shown in Q, however, that the features in Fig. 2 
(most notably the v spectral split) may be smeared out in 
more realistic ("multi-angle" Q) numerical experiments 
accounting for different (unaveraged) trajectory angles. 

Here wc study the transition from single- to multi- 
angle simulations in a (hypothctically) more favorable 
scenario, where (E x ) is decreased from 24 MeV to 18 
MeV — everything else being unchanged. This choice, due 
to a "fatter" low-energy tail, leads to a slightly more 
prominent v split at higher- E c , as evident in the single- 
angle results of Fig. 5 (to be compared with Fig. 2). 

Concerning the reduction to four modes, we find that 
the case in Fig. 5 is well described by E c ~ 8.2 MeV, 
~E C ~ 6.1 MeV, and: J, ~ 0.131, J h ~ 0.031, Ji ~ 0.012, 
J/j ~ 0.029. The corresponding frequencies are: u>i ~ 
1.16, u) h ~ 0.97, ZJi ~ 1.36, and uJ h ~ 0.66 (in km" 1 ). 

Concerning the transition from single- to multi-angle 
calculation, however, we find again that the latter tend 
to smear out the low-energy spectral features. Figure 6 
(multi-angle) shows, in comparison with Fig. 2 (single- 



angle) that: (1) the v split is broadened, and (2) the V one 
is largely suppressed, and survives as a slight low-energy 
"shoulder." At present, this (unfavorable) numerical ob- 
servation [l7T ] remains analytically unexplained. One can- 
not exclude, however, that the v spectral split feature 
may be less suppressed in other multi-angle cases, whose 
investigation is left to future work. 

V. SUMMARY 

We have studied some low-energy features of the V en- 
ergy spectrum, focusing on the spectral split phenomenon 
emerging (in inverted hierarchy) in single-angle simula- 
tions. We have related its origin to nonadiabatic aspects 
of the v and V evolution, which can be simply modeled 
through an alignment ansatz with four energy modes. 
However, the V split feature appears to be fragile when 
passing to multi-angle simulations. Further studies are 
required to deepen its analytical understanding, as well 
as the conditions for its observability. 
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FIG. 6: As in Fig. 5, but in multi-angle simulation. 
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